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Abstract. In the present paper we are concerned with the Novikov-Veselov equation at 
negative energy, i.e. with the (2 + 1) -dimensional analog of the KdV equation integrable by the 
method of inverse scattering for the two-dimensional Schrodinger equation at negative energy. 
We show that the solution of the Cauchy problem for this equation with non-singular scattering 
data behaves asymptotically as "^"^^ in the uniform norm at large times t. We also prove that 
this asymptotics is optimal. 

o. 

^_, : 1 Introduction 

^J ! In the present paper we consider the Novikov-Veselov equation 

^ : dtv = ARe{Ad^,v + d^{vw) - Ed,w), 

pT"; d2W = -3d^v, v = v, EgR, E < 0, (1.1) 

■^ . V = v{x,t), w = w{x,t), X = (xi,X2) G ffi , t £ M, 

^_> , where 

C^ ■ „ d ^ 1 f d . d \ ^ \ f d . d 



X 



dt' ^ 2 \dxi 8x2) ' ^ 2 \dxi 8x2 

We will say that {v,w) is a rapidly decaying solution of (jl.ip if 
CN ■ 

^: •v,w€C{R'^ xR), v{-,t) eC^{R^), (1.2a) 

l/^ ! a(t) 

- ^div{x,t)\ ^ ^^ , \j\ ^ 3, for some e > 0, w{x,t) -^ 0,\x\ -^ 00, (1.2b) 

(1 + |x|) 

r~>^ . • (v^w) satisfies (jl.l|) . (1.2c) 

o: 

Note that if v{x,t) = v{xi,t), w{x,t) = w{xi,t), then p.ip is reduced to the classic 
KdV equation. In addition, (jl.ip is integrable via the inverse scattering method for the two- 
dimensional Schrodinger equation 



Lip = Eip, L = -A + v{x), x = {xi,X2), E = Efixed- (1-3) 

In this connection, it was shown (see [Mj, |NV1] . [NV2] ) that for the Schrodinger operator L 
from (|1.3p there exist appropriate operators A, B (Manakov L-A-B triple) such that (jl.ip is 
equivalent to 

^-^^^ = [L-E,A] + B[L-E), 

where [•, •] is the commutator. 

Note that both Kadomtsev-Petviashvili equations can be obtained from (jl.ip by considering 
an appropriate limit E — )• itcxD (see |ZS] . [G2j ). 

We will consider the Cauchy problem for equation (jl.ip with the initial data 

v{x,{)) = vo{x), w{x,Q) = wq{x). (1-4) 
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We will assume that the function ^0(2;) satisfies the following conditions 

• 7;o = vo, (l-5a) 

• £^ = 0, where E is the set of zeros of the Predholm determinant A of (|2.7p 
for equation (|2.6p with v(x) = Vo{x), 

• f G 5(M ), where S denotes the Schwartz class. (1.5c) 



(1.5b) 



As for the function W(){x), which plays an auxiliary role, we will assume that it is a continuous 
function decaying at infinity and determined using dzWo{x) = —3dzVo{x) from (jl.ip . 

Condition (jl.Sbp is equivalent to non-singularity of scattering data for ^0(2;). Conditions 
(jl.Sp define the class of initial values for which the direct and inverse scattering equations ()2.4p , 
(|2.9p - (|2.1ip . with time dynamics given by (|2.15p . are everywhere solvable and the corresponding 
solution V of (jl.ip belongs to C°°(M^, M). We will call such solution {v{x, t),w{x, t)), constructed 
from {v(){x),wq{x)) via the inverse scattering method, an "inverse scattering solution" of (jl.ip . 

The main result of this paper consists in the following: we show that for the "inverse 
scattering solution" v{x,t) of (jl.ip . (jl.4p . where E < and v{x,0) = vo{x) satisfies (jl.Sp . the 
following estimate holds 

const(wo) ln(3 + Itl) _ , , _9 , , 

\v{x,t)\ ^ ^ "^ ^3 ' '^ t G M, uniformly on x G M^ (1.6) 

We show that this estimate is optimal in the sense that for some initial values v{x, 0) and for 
some lines x = ut, cj G S^, the exact asymptotics of v{x,t) along these lines is , "|'"^3/4 as 
1^1—7-00 (where the constant is nonzero). Note that de facto the "inverse scattering solution" 
is the rapidly decaying solution in the sense of (jl.2p . 

This work is a continuation of the studies on the large time asymptotic behavior of the 
solution of the Cauchy problem for the Novikov-Veselov equation started in [KNlj for the case 
of positive energy E. It was shown in [KNl] that if the initial data {vo(x), wo{x)) satisfy the 
following conditions: 

• {vq{x), wo{x)) are sufficiently regular and decaying at |x| — ;■ 00, 

• vo{x) is transparent for l\l.'d\] at E = Efi^ed > 0, i.e. its scattering amplitude / is 
identically zero at fixed energy, 

• the additional "scattering data" b for vq{x) is non-singular, 

then the corresponding solution of (jl.ip . (jl.4p can be estimated as 

, , ,, const • ln(3 + Itl) ™ ^ , ^9 

\v{x,t)\ ^ — — -^, t G M uniformly on a; G R . 

' \ \ 

This estimate implies, in particular, that there are no localized soliton-type traveling waves 
in the asymptotics of (jl.ip with the "transparent" at E = E fixed > Cauchy data from the 
aforementioned class, in contrast with the large time asymptotics for solutions of the KdV 
equation with reflectionless initial data. 

It was shown in [N2j that all soliton-type (traveling wave) solutions of (jl.ip with E > Q 
must have a zero scattering amplitude at fixed energy; in addition it was proved in [N2] that 
for the equation (jl.ip with E > Q no exponentially localized soliton-type solutions exist (even 
if the scattering data are allowed to have singularities). However, in |G1] . [G2j a family of 
algebraically localized solitons (traveling waves) was constructed de facto (see also |KN2j ). We 



note that for the case E < 0, though the absence of exponentially-locahzed sohtons has been 
proved (see jKN3j ). the existence of bounded algebraically localized solitons is still an open 
question. 

Note that studies on the large time asymptotics for solutions of the Cauchy problem for the 
Kadomtsev-Petviashvili equations were fulfilled in |MST| . |HNS| . p?] . 

The proofs provided in the present paper are based on the scheme developed in |KNlj . the 
stationary phase method (see [Fej ) . These proofs include, in particular, an analysis of some 
cubic algebraic equation depending on a complex parameter. 

This work was fulfilled in the framework of research carried out under the supervision of 
R.G. Novikov. 

2 Inverse "scattering" transform for the two-dimensional Schro- 
dinger equation at a fixed negative energy 

In this section we give a brief description of the inverse "scattering" transform for the two- 
dimensional Schrodinger equation ()1.3p at a fixed negative energy E (see |GN| . |N1| . |G2j ). 

First of all, we note that by scaling transform we can reduce the scattering problem with an 
arbitrary fixed negative energy to the case when E = —1. Therefore, in our further reasoning 
we will assume that E = —1. 

Let us consider potentials v{x) for the problem ()1.3p satisfying the following conditions 

^ = ^' 1^(^)1 ^ (l + |l|)2+e ' ^e^^'' (2-1) 

for some fixed q and e > 0. Then it is known that for A G C\(0 U £), where 

£ is the set of zeros of the modified Fredholm determinant A 

„ , (2.2) 

for equation (j2.6p . 



there exists a unique continuous solution tl^{z, A) of (jl.Sp with the following asymptotics 

V'(z,A) =e-^(^^"+^A)/i(z,A), /x(z,A) = l + o(l), \z\ ^ oo. (2.3) 

Here the notation z = xi + ix2 is used. These solutions are known as the Faddeev solutions for 
the Schrodinger equation ()1.3p , E = —1, see for example |Fa) , |Nlj . 
The function fi{z, A) satisfies the following integral equation 

Kz, ^) = ^ + jj 9{^ - ^' ^XOKC A)dReC(iImC (2.4) 

.(.,A) = - f J-V //exp(V2(a + C.)) 

""^ ' \2ttJ J J CC + i(AC + C/A) ^ ' 

Cec 

where z G C, A G C\0. 

In terms of m^z, A) = (1 + |2;|)~'^^^''^/i(z. A) equation (|2.4p takes the form 

m(z. A) = (1 + |z|)-(2+.)/2 + ||(i + |^|^_(2+.)/2^(^ _ ^^ A)__ZiO__^(^^ A)dReCdImC, 

Cgc 

(2.6) 



where z G C, A G C\0. In addition, A{-,-,X) £ L^{C x C), |Tr^2(A)| < oo, where A{z,CA) 
is the Schwartz kernel of the integral operator A{X) of the integral equation (|2.6p . Thus, the 
modified Fredholm determinant for ()2.6p can be defined by means of the formula: 



In A(A) = TY(ln(/ - A{X)) + A{X)) (2.7) 

(see (_GK ] for more precise sense of such definition) . 

Taking the subsequent members in the asymptotic expansion (j2.3p for -0(2;, A), we obtain 
(see [Nlj): 

ip{z, A) = exp f -- (^Az + I M ll- 27rsgn(l - AA) x 

\z\ ^oo, XeC\{£LlO). 

The functions a(A), 6(A) from (j2.8p are called the "scattering" data for the problem ()1.3p . 
(I2T]) with E = -1. 

The function fj,{z, A), defined by (j2.4p . satisfies the following properties: 

/^(z. A) is a continuous function of A on C\(0 U £); (2-9) 

dn{z,X) 



dX 



r{z,X)n{z,X), (2.10a) 



r(z, A) = r(A)exp Q (^(^A - i) z - (a - ^) z)) , (2.10b) 

vrsgnfl — AA) , ,, , , , 

r(A) = ^ \ -b{X) (2.10c) 

A 



for AGC\(OUf); 



^ — > 1, as A — > oo, A -^ 0. (2-11) 

The function 6(A) possesses the following symmetries: 

1\ /I 



6{-=)=6(A), 6( = )=6(A), AgC\0. (2.12) 



In addition, the following theorem is valid: 
Theorem 2.1 (see [GN], [NT], [G2]). 



(i) Let V satisfy /ll.5\) . Then the scattering data 6(A) for the potential v{z) satisfy properties 
^KM "-nd b G S{D^), where L>_ = {A G C: |A| > 1}, Z)_ = L>_ U dD^ and S denotes 
the Schwartz class. 

(a) Let b be a function on C, such that 6 G 5(l)_) and the symmetry properties Ii2.12\) hold. 
Then the equations of inverse scattering k2.1b]) - l2.18\) are uniquely solvable and the cor- 
responding potential v{z) satisfies the following properties: v G C°°(C), v = v, \v{z)\ — )• O 
when \z\ — )• oo. 



Let us denote by T the unit circle on the complex plane: 

r = {AGC: |A| =1}. (2.13) 

Then, in addition, it is known that under the assumptions of item (i) of Theorem 12.11 the 
function h is continuous on C and its derivative d\b{\) is bounded on C, though discontinuous, 
in general, on T. 

Finally, if {v{z,t),w{z,t)) is a solution of equation (|1.1|) with E = —1 in the sense of (|1.2|) . 
then the dynamics of the "scattering" data is described by the following equations (see [GNJ ) 

a(A,t) = a(A,0), (2.14) 

6(A, t) = exp I ("a^ + 1 - P - 1") t\ 6(A, 0). (2.15) 

The reconstruction of the potential v{z, t) from these "scattering" data is based on the 
following scheme. 

1. Function /x(z, A,t) is constructed as the solution of the following integral equation 

M^,A,t) = 1 -l.ll^r{z,C,t)WxJ)^^^^- (2.16) 

2. Expanding /i(z. A, t) as A — )• oo, 

M-,A,t) = l + ^^^i^+o(^), (2.17) 

we define v{z,t) as 

v{z,t) = -2d,iJ.i{z,t). (2.18) 

3. It can be shown that 

Lip = El/; 

where 



L = -4d,d2 + v{z,t), v{z,t) = v{z,t), E = -l 
ijiz, A, t) = e-^(^^+^/^V(^, A, t), A e C, zeC, te 



3 Estimate for the linearized case 

Consider 



I{t, z) = II /(C) exp(5(C, z, t))dReCdlmC, 

c 

J{t, z)=^ll ^/(C) exp(5(C, z, t))dReCdImC, 



where z G C, t G M, /(C) G L^C), S is defined by 



1 //. 1\ /. 1\ \ /., 1 ^, 1 

=3 



S{X,z,t) = ^([\-j.jz-^\-j-jzj+t^\^ + j^-X^-^]. (3.2) 



(3.1) 



We will also assume that /((") satisfies the following conditions 

/GC~p+), fGC^iD.), (3.3a) 

where 

i?+ = {AGC: 0< |AK 1}, D^ = {\eC:\X\>l}, (3.4) 

and Dj^ = D+UT, D^ = D^UT with T defined by (i2J3]l . 
Note that if v(z, t) = /(t, z), w(z, t) = J{t, z), where 

(ICI^ + ICr^)/(C) G ^^(C) as a function of C, 
and, in addition, 

4„n 



/(C) = /(-C) and/or f{C) = -\C\-f^-. 
then f , ty satisfy the linearized Novikov-Veselov equation (jl.ip with i? = — 1. Besides, if 

/(C) = ""'^^"if KC), (3.5) 

where 6(C) is the scattering data for the initial functions {vq{z),wo{z)) of the Cauchy problem 
(jl.ip . ()1.4p . then the integrals /(t, z), J{t, z) of ()3.ip represent the approximation of the solution 
{v{z,t),w{z,t)) under the assumption that || v ||<C 1. 

The goal of this section is to give, in particular, a uniform estimate of the large-time behavior 
of the integral I(t, z) of (|3.ip . 

For this purpose we introduce parameter tt = | and write the integral / in the following 
form 

I{t,u)= fff{C)exp{tS{u,C))dReCdlmC, (3.6) 

c 

where 

5(».C) = ^(c-i)^-(c-4)») + (c»-? + i-^). (3.7) 

We will start by studying the properties of the stationary points of the function S{u, Q with 
respect to C- These points satisfy the equation 






S'c = 7,-l^ + K'-iz = 0- (3. 



The degenerate stationary points obey additionally the equation 



u , 12 
We denote C = C^ and 



^CC = 73+6C + 75=0- (3-9) 



QK0 = 7T-;T7 + 3C c2- 



For each ^, a root of the function Q{u,S,), there are two corresponding stationary points of 

5(n,C),C = ±\/?. 



The function SUuX) can be represented in the fohowing form 

s'^{u,o = ^(c' - ciinm' - diume - cHu)). (s.io) 

We win also use hereafter the fohowing notations: 

U = {u = -6{2e-''^ + e^'^P) : 99 G [0, 2^)} 
and 

U = {n = re^^: V' = Arg(-6(2e-^^ + e^^^)), O^r^ |6(2e^'^ + e'^*'^)], 99 G [0,27r)}, 
the domain hmited by the curve lA. 
Lemma 3.1 (see [KNT]). 

1. Ifu = -18e^, k = 0,1,2, then 

Co{u) = Ci{u) = C2{u) = e^ (3.11) 

and S{u, C) has two degenerate stationary points, corresponding to a third-order root of 
the junction Q{u,t,), ii = e 3 . 

2. IfueU{ i.e. u = -6(2e-*'^ + e^*^) ) and u ^ -18e^, A: = 0,1, 2, then 

Co(n) = Ci{u) = e-'^l\ C2(n) = e'^ . (3.12) 

Thus S{u, C) has two degenerate stationary points, corresponding to a second-order root of 
the function Q{u,^), ^1 = e"**^, and two non-degenerate stationary points corresponding 
to a first- order root, ^2 = e'^'^'^. 

3. Ifue intU = U\aU, then 

Ci{u) = e"*"^' for some real ipi, and Cii^) 7^ Cjiu) for i ^ j. (3.13) 

In this case the stationary points of S{u,Q are non- degenerate and correspond to the roots 
of the function Q{u,^) with absolute values equal to 1. 

4. Ifue C\V, then 

Co(tx) = (l+a;)e-*^/2^ Ci(n) = e*^ C2{u) = {1 + u^ye-'^P/^ (3.14) 

for certain real values 93 and uj > 0. 

In this case the stationary points of the function S{u, C) are non- degenerate, and cor- 
respond to the roots of the function Q{u,^) that can be expressed as ■^o = (1 + T)e~^'^, 
6 = e'^'^P, ^2 = (1 + ry^e-'^P, (1 + r) = (1 + uf. 

Formula (j3.10p and Lemma 13.11 give a complete description of the stationary points of the 
function S{u, (). 

In order to estimate the large-time behavior of the integral having the form 

I{t, u, A) = ff /(C, A) exp{tS{u, C))dReCdlmC (3.15) 

c 

(where S{u, Q is an imaginary-valued function) uniformly on u,X £ C, in the present and the 
following sections we will use the following general scheme. 



1 . Consider D^ , the union of disks with a radius of e and centers in singular points of function 
f{C-,^) and stationary points of S{u,C,) with respect to C,. 

2. Represent I{t,u,X) as the sum of integrals over Z^^ and C\L'e: 

I{t, u, A) = lint + lext, where 

I^nt = If /(C, A) exp(t5(n, C))dReCdImC, 

ii (3.16) 

hxt = II fiC, A) exp(tS(n, C))dReCdImC. 

3. Find an estimate of the form 

\Iint\ = O {e'') , as e^O (q ^ 1) 
uniformly on u, A, t. 

4. Integrate lext by parts using Stokes formula 

^"*" 2zt J S'inX) ^ 

1 f (/.(C, A) - /, (C, A)) exp(t^(n, Q) ,^-_ 1 ^^ /c(C, A) eMtSju, Q) ^^^^^^^^^ 



2it 7 «(n,C) tyy 5'(n,C) 

'' ^^ dReC(iImC = --(Ii+/2 + /3--f4), (3.17) 



tJJ (S'^iuX)? ' ' i 

C\D, 

where /±(C, A) = lim /(C(l T S), A), C e ^ and T is defined by pTTB]) . 

5. For each Jj find an estimate of the form 

\Ii\ = O I —g j , as e — ^ 0. 

6. Set e = -^ — \j.\\k ' '^^^sre k{a + /3) = 1, which yields the overall estimate 



1 

(l + |t| 



\I{t, u, X)\ = O { —^- , as \t\ — )• oo. 

V(i + |i|)"+V 

Using this scheme we obtain, in particular, the following result 
Lemma 3.2. Let a function f satisfy assumptions l[il.&]) and, additionally, 

/|t = 0, r = {AGC:|A| = l}. (3.18) 

Then the integral I{t,u) of i3.6\) can be estimated 

uniformly on u £ C 

Note that condition ()3.18p is satisfied if / has the special form (|3.5p . A detailed proof of 
Lemma 13.21 is given in Section [6l 



4 Estimate for the non— linearized case 

In this section we prove estimate (jl.6p for the solution v{x,t) of the Cauchy problem for the 
Novikov-Veselov equation at negative energy with the initial data v(x,0) satisfying properties 

dESD. 

We proceed from the formulas ()2.17p . ()2.18p for the potential v{z, t) and the integral equation 
(|2J6]1 for ^(z,A,t). 
We write (|2J6]l as 

fi(z,X,t) = l + {A,^tf^)(z,X,t), (4.1) 

where 



c 



2^ 

and S{uX) is defined by (j3.7p . u= —. 



Equation ()4.ip can be also written in the form 

H{z, A, t) = 1 + A,,t • 1 + {Altl^){z, A, t). (4.2) 



According to the theory of the generalized analytic functions (see [V]), equations ()4.ip . (|4.2p 
have a unique bounded solution for all z, t. This solution can be written as 

f,iz,X,t) = iI-Al,y\l + A,,t.l). (4.3) 

Equation ()4.3p implies the following formal asymptotic expansion 

l,{z, A, t) = (/ + A^^, + 4, + ...)(! + A,,i • 1). (4.4) 

We also introduce functions u{z,X,t) = dzfJ-{z,X,t) and r]{z,X,t) = dzn{z,X,t). In terms of 
these functions the potential v{z,t) is obtained by the formula 

v{z,t) = -2u_i{z,t), (4.5) 

where z^_i(z, t) is defined by expanding ^{z, X,t) as |A| — )• oo: 

The pair of function i^(z, A, t), ri{z, X, t) satisfy the following system of differential equations: 
dv{z^X,t) 



dX ^ 
dr]{z,X,t) 



dzr{z, X, t)fi{z, A, t) + r{z, A, t)r]{z, A, t) , 



(4.6) 



dX 
Equations ()4.6p can also be written in the integral form 



dzr{z, A, t)/x(z, A, t) + r{z, A, t)z^(2;, A, t) . 



u{z, A, t) = {Bz^t^j){z., A, t) + {Az^tri){z, X, t), 
r]{z, A, t) = {Bz,tf^){z, A, t) + {Az,tT^){z, A, t), 



(4.7) 



where operator B^ t is defined 



{Bz,tf){X) = d^Hd,riz,X,t)fiX)) = -lyyMlMZKldReCdlmC. (4.8) 



Thus for the function iy{z,X,t) we obtain equation 

or the following formal asymptotic expansion 

z. = (/ + Al, + AX, + . . .)((i?.,t + A,^tB,,t){I + A% + 4, + ...)(! + A,,t • 1)). (4.9) 

We will write this formula in the form v = B^^t ■ 1 + Az^tBz,t ■ 1 + Rz,ti^)- 
Lemma 4.1. Let f{X,z,t) be an arbitrary testing function such that 



[i + \t\Y' l^^-^l ^ (1 + |t| 



I/KttTW' I^^-^I ^TTtW VAGC,zGC,tG 



with some positive constant cj independent of X, z, t and some 6^0. Then: 
1. The following estimates hold for B^^t " /■' 

iBz,t ■ /)(A) = ^^ + a (^) for X ^ oo, 



where 



and 



Pi{z,t) 



^ 11 dzr{z, C, t)/(C, z, t)dReCdlmC, 



in addition, 



'"'•^•''KlTTM^^ <"»' 



(B.,./)(A)i< w+iy ,„,,,^. „,,, 



where T is defined by Ii2.1°:l\) . and 



l(^^--/)(^)K(T^q|W ^^^^- (^-12) 



S. The following estimates hold for Az^t • B^^t • f ■ 



{Az,t ■ Bz,t • /)(A) = ^i^ + o (^ ) /or A ^ oo, 



where 



aiiz,t) = -\ ff dReCdlmCr{z,C,t) If ""^^^'^'^^ fiv,z,t)dRer]dlmrj, (4.13) 



and 



in addition, 



l«i(^'^)l^ (l +",1)3/4+. ; (4-14) 



\{Az,fBz,ff){X)\^ (l+'lt'oW^ ^^^'^^ (^-^^^ 
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3. The following estimates for A^^ ■ f hold: 



where 



Kr/)(A) = ^^ + o(^) /orA^oo, 



1 



7„(z,t) = -JJ rizX,t){A^,y ■ fXOdReCdimC 
c 
and 

\^n{z,t)\ ^ }^'lll.,.. , (4.16) 

where [s] denotes the smallest integer following s. In addition, 

|(Ag,-/)(A)K ^^'jy^r., VAGC. (4.17) 

4- The following estimate holds for Rz,t- 

Rz,tW = %^ + o (^) for A ^ oo, (4.18) 

and 

A detailed proof of Lemma l4. II is given in Section [6l 

Prom formulas ()4.5p . ()4.9p and Lemma |4. II follows immediately the following theorem. 

Theorem 4.1. Lei v{x,t) be the "inverse scattering solution" of the Cauchy problem for the 
Novikov-Veselov equation lil.l]) with E = —1 and the initial data w(a;,0) = vq{x) satisfying 
nrm . Then 

, , ,, const (wn) ln(3 + lil) ^o ^ 

\v{x,t)\^ (1 + 1^1)3/4 ' xGM2,tGM. 

5 Optimality of estimates (11.61) and ( 13.191) 



In this section we show that estimates ()1.6p and ()3.19p are optimal in the following sense: 
there exists such a line z = iit that along this line I{z, t) from (|3.19p behaves asymptotically 
as ^-1 ?°HN3/4 with come nonzero constant; there exist such initial data satisfying (|1.5p that the 

corresponding solution v{z, t) behaves asymptotically as ,/^°"'?3/4 as |t| — )• oo, where the constant 
is nonzero. 

5.1 Optimality of the estimate for the linearized case 

Let us consider the integral 



I{t, u) = II /(C) exp(t5(n, C))dReCdImC, (5.1) 



where 

vr|l-CCI 

2|Cl 



/(c) = '.,.,r' Kc), (5.2) 
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with 6(C) satisfying ()3.3p . and S{u, Q is defined by ()3.7p . for u = n = — 18. As shown in Lemma 
13. II for this value of parameter u the phase S{u, C) = S{Q has two degenerate stationary points 
(^ = ±1 of the third order. 

To calculate the exact asymptotic behavior of I(t, u) we will use the classic stationary 
method as described in [E]- First of all, we note that /(C) is continuous, but not continuously 
differentiable on C. Thus we will consider separately the integrals 

I+(t) = II /(C) exp{tS{0)dReCdImC, I^{t) = I I /(C) exp(t5(C))dReCdImC, 

D+ -D- 



where -D_|_ and D- are defined in 

Let us introduce the partition of unity V'i(C) + ''Ao(C) + ^^-1(0 = li such that ^ Vj ^ Ij 
V'j G C°°(C), '(A±i(C) = 1 in some neighborhood of C = il, respectively, and V'±i(C) = 
everywhere outside some neighborhood of C = il respectively. Then it is known (see [^) that 

I+{t) = II f{C)MC)expitS{C))dReCdImC+ll f{C)^P-i{C)ew{tS{C))dReCdlmC+0 (^^ 

= I+{t) + I^{t) + o(^/\ as \t 



D+ D+ 

oo. 



First we will estimate I^{t). We note that for the phase ^(C) the following representation 
is valid 

S{C) = P{C) - P{C), 

where -P(C) is a holomorphic function defined by 

P(C) = -9C-^ + C' + ;^ + l6; 

in addition, P^{C) = 5^(0 = ^(C - 1)'(C + !)'• 

Note that P(C) can be written in the form P(C) = /'(C) (C ~ 1)^) where p{C) = a3 ^.nd 
lim /9(C) 7^ 0. For the function p{C,) the expression {piOY'^ can be uniquely defined in some 
neighborhood of C = 1- Further, we define the transformation (^ ^ rj: 

Since we have that 



drj 



= ^/0, (5.3) 

C=i 

the inverse transformation C = '^{v) is defined in some small neighborhood of r/ = 0. In terms 
of the new variable rj the phase can be represented 



Now if we denote x = Kerj, y = Imr/, the integral I^{t) becomes 
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where f = f o ip^ tfji = t/ji o ip and A-j. = {{x,y) € M? : x < 0}, i.e. A+ is the half-plane 
containing the image of -D+ n Bgr^l) under the transformation x + iy = ip~^{(). 
The integral /+(t) can be written in the form 

+00 



^+(*) = / dcexp(3itc) / f{x + iy)%l^i{x + iy)\dr,(p{x + iy)\ dujs, 

-00 7cnA+ 

where dus is the Gelfand-Leray differential form, defined as 

dS A dios = dx Ady (5.4) 

and in the particular case under study equal to 

-{x^ — 3xy'^)dx + (3x^y — y^)dy 



dujs 



{x? + y 



2^3 



7c is an oriented contour consisting of points of the set {S{x, y) = c} with the orientation chosen 
so that (I53D holds. 

As ijjiix + iy) is equal to zero outside some Br(0), a disk of radius R centered in the origin, 
then there exists such c* > that the set {S{x,y) = c} lies outside -B_r(0) for any c < — c#, 
c> c^,. Thus the integral I^{t) can be written 



I^{t) = / dcexp{3itc) / f{x + iy)ipi{x + iy)\dri(p{x + iy)\ dojs- 

-c, 7cnA+ 

Performing the change of variables 

{X I—)- c^'^x, y I—)- c^i^y for c > 0, 
X I-)- {-cY^'^x, y H> (-c)^/^y for c < 



yields 



where 



c. 



/■ dcexp(3itc) „ . . /" dcexp(3itc) 







F+(c) = / J='{c,x,y)dujs, ^-(c) = / T{-c,x,y)dws, (5.5) 

7+, 7_ are oriented sets consisting of points of the sets {S{x,y) = 1}, {S{x,y) = —1} corre- 
spondingly with orientation chosen so that ()5.4p holds and 

He, X, y) = ~f{c^l\x + iy))V;i(ci/^x + iy))|a,^(ci/^(x + iy))\\ 

For any fixed positive c the integrals in ()5.5p converge because the set {5(x, y) = 1} is separated 
from zero and, consequently, the denominator does not vanish, and because Vi is a function 
with a bounded support and thus the domains of integration in (|5.5|) are, in fact, bounded. 
Besides, since A+ is a conic set, the functions -F+(c) and F^{c) can be expressed as follows 



F+{c) = j ^{c,x,y)dus, F_{c)= LF{-c,x,y)dujs, (5.6) 

7+nA+ 7_nA+ 
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In some neighborhood Uq of c^'^(x + iy) = containing the support of -01 the function 
f{c^'^{x + iy)) can be represented as 

fic'/\x + iy)) = /(I) + 6-i/4[5c/i,+ (l)(x + iy) + 5^-/d+(1)(x - iy)]c'/'' + g{c'/\x + iy)), 

where dc^fo+iX) — 1™ ^c/' ^c-^£>+(l) = 1™ "^j/ and, in addition, (7 is a function that can be 

estimated 

{c^l\x + iy))\ ^ K\c^l\x + iy)|i+°, ci/4(x + iy) G ZYq (5.7) 



|5 



def 



with some constants a > 0, i^ > 0. Using ()5.2p . we note that /(I) = 0, 9^/d+(1) = <9a/d+(1) 
/^ (1) and thus 

f{c^l\x + iy)) = Q~^'^f'j,^{l)xc^l^ + g{c^'\x + iy)), c^/\x + iy) eU,. 
Taking into account (|5.3p we obtain 

F{c,x,y)=^f'ij^{l)x^'^ + ~g{c^'\x + iy)), c^'\x + iy) eU^, (5.8) 

where 7 = 6~^'^ and 5 is a function satisfying an estimate similar to ()5.7p . 
It fohows then that the functions F±{c) behave asymptotically as 

i^±(c) = 7/D+(l)^A+(±c)'/' + ii(c), whenc^O, 

where 

7+nA+ 7_nA+ 

and R{c) denotes the remainder. 

The integrals J^ converge because the set {S{x,y) = ±1} represents a combination of 
curves which do not pass through zero and converge either to the lines \y\ = \x\ or to the 
coordinate axes with velocities \y\ = tA^ and |x| = |4^ correspondingly. 

The remainder R{c) behaves asymptotically as o (c^' ^) because we can estimate 

|/?(c)| ^ i^ci/^(i+") / \x + iy\^+'^\dws\. 

(7+U7-)nA+ 

and the integral converges due to the properties of the set {S{x,y) = ±1} explained above. 
Thus I^{t) behaves asymptotically as (see [El Chapter III, §1]) 



n(t) = ^J'o.mr [^ 



r - '"' 



Ja. exp ^ ^ + Ja , exp ( — 



^3/4+^(^^3/4 J 



where F is the Gamma function. 

Let us perform the same procedure for I^{t). In this case we will define P{C) as -P(C) 
-9C - 7 + C^ + 7^ - 16 and obtain 



C ' ^ ' C^ 



P{C) = p{C){C + i)\ p(C) = ^^^^. 
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Since />(— 1) = —6, we will define the following transformation C — )• r/ in the neighborhood of 
( = —1: rj = {—piOYiC + 1) (note that ^ = v^). Then 5(C) can be represented 

5(C) = —rj^ + f]^ and integral I^{t) becomes 

A_ 

where A„ = {{x,y) €M?: x > 0} and the functions /, "0-1) ^ are defined similarly to the case 
of I^{t). The integral I^{t) can also be written 



/_^(t)= / dcexp(-3itc) / f{x + iy)%lji{x + iy)\d.^ip{x + iy)\ dujs, 

— oo 7cnA_ 

where 7c and (ia;5 are the same as for the case of I^{t). Performing further the same procedure 
as for the case of I^{t) and taking into account that J^ = —J^ , J^ ~ ~'^a ' ^^ obtain 



the following asymptotic expansion for I^{t): 



-^A+expl-— j+J^^expj^- 



^+'^f^)- 



Considering the case of /_ (t) we note that in order to get an asymptotic representation for 
I^{t) and lZ{t) we need to replace Z)+ — t- Z)_, A^- — t- A_, A_ — t- A+ in the formulas for I^{t) 
and /+ (t) correspondingly. Taking into account that f'^ (1) = — /^_ (1), f'f) ( — 1) = — /£>_ ( — !)> 
we obtain 

where 

--|^^(!)(/^.«(^t-(")--^.-(-"))- 

/d,(-1)(4,«p(-^)+'^a,«p(^))). (5.9) 

Thus we have shown that the linear approximation of the solution v{z,t) of (jl.ip . when 
z = — 18t, behaves asymptotically as -^ when t — ;■ oo. 

Note that on the set 7+U7_ the differential form dujs is positive. Thus J^ are some negative 
constants, and expressions J^ exp (^) + J^ exp(-^), J^ exp (-^) + J^ exp (^) 
do not vanish. On the other hand, from (j2.12p and (j5.2p it follows that f'j^ (l)//j-, (—1) = 

6(1)/5(1). Thus, in the general case the constant C from ()5.9p is nonzero. We have proved the 
optimality of the estimate p.6p in the linear approximation. 

5.2 Optimality of the estimate for the non-hnear case 

Now we show that for certain initial values v{z,0) the corresponding solution v{z,t) of (jl.ip 
behaves asymptotically as -^ along the line z = —18t for some c / 0. Let us show that the 
integral (|4.13p with f =1 and z = -18t behaves as ^yf^. 
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We can represent qi(— 18t,i) in the form ()5.ip with f{C,t) = '^(C;0/'(C)0) where 

P(C, = -i // ^^^^^^^f^^dRevdReC 

c 



In other terms, 

/(C) = /(C,t) = ^1^^(^6(00.^ Ni-cci 



6(C)exp(t5(C)) 



2|CP 
We proceed following the scheme of estimate for I{t) until formula ()5.7p . Then we represent 

V;i(ci/^(x + iy)m,ip{c^/\x + iy))|2 = ^ + fcci/^(x + iy) + /i(c^/^(x + iy)), 

V6 

where A: is some coefficient and h{c^'^{x + iy)) satisfies an estimate of type ()5.7p . Consequently, 
for J-'{c, X, y) we can write 

T{c,x,y) = f{l){6-'/^+kc'/\x+iy))+6-'/'[d^fD4^){x+iy)+d^fD4l){x-iy)]c'^ 
where g{c^''^{x + iy)) satisfies an estimate of type (j5.7p . This allows us to obtain in the end 
r i«, ,^ /i(/(±l)) f, , const^ , /2(/c(±l),/c(^l)) , r 1 A ,^ ,,in^ 

where /i(/(ibl)) is a linear combination of the limit values of / as C tends to 1 and —1 from 
inside and outside of the unit circle, and l2{f({±l), f?{±l)) is a linear combination of the limit 
values of /^ and /^ as C tends to 1 and —1 from inside and outside of the unit circle. 

Now let us consider the potential vg corresponding to the scattering data 9b, where ^ G M 
is some small parameter. In a way similar to which ()5.10p was obtained it can be shown that 

|/(±l,t)| < 1^, |/c(±l,t)| ^ C29^ |/f(±l,t)| ^ €20^ 

for sufficiently large values oft, where ci, C2, C3 are some constants independent of i and 6. When 
0—7-0 and t — >• 00, the linear approximation of vg behaves as O i -374 1 , while the expression 

ai(— 18i, t) behaves as O i -^ J (it can be shown that the member o i -Jj^ J in ()5.10p depends 
quadratically on 9). 

Finally, from ()4.9p and Lemma l4. II it follows that for 9 small enough and for z = — 18i, 

ve{z,t) = ^ + o^-^j, t^oo, 

where Cq is some nonzero constant. Thus we have shown that the estimate ()1.6p is optimal. 

6 Proofs of Lemma 13.21 and Lemma 14.11 

Proof of Lemma \3.S\ The proof follows the scheme described in Section [3] and is carried out 
separately for four cases depending on the values of the parameter u. In all the reasonings that 
follow we denote by Dg the union of disks with the radius e centered in the stationary points 
of S{u, Q) and we denote by T the unit circle on the complex plane: 

r = {AGC: |A| =1}; (6.1) 
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in addition, const will denote an independent constant and const (/) will denote a constant 
depending only on function /. 

Case 1. uGV 

In this case all the stationary points belong to T and due to assumptions (|3.3p and (|3.18p 
of Lemma 13.21 we can estimate 



■2) 



1/(01 ^ const(/)e for C G D^ 
Now we estimate the integral lint (as in ()3.16p ) as follows 



'iratl 



Jl /(C) eMtS{u, C))dReCdlmC 



Oh 



^ const (/) e dReCdImC ^ const (/)e^. 



D, 



The estimate for lext (as in (j3.16p ) is proved as follows. 
We note that the function SUu, C) can be estimated as 



3 



\S[{u,0\ > const -5^ for C G C\D,„ and 



|5^(n,C)| ^ const T— J for C G dDp, e ^^ p ^ Eq. 



Similarly, we can estimate 






iS'Au,c)f 



^ const — J- for C G C\L'e(,, and 



|C|4 
^ const — ^ for ( £ dDp, e ^ p ^ Eq. 



Thus we obtain the following estimate for Ii from (j3.17p 



(6.3) 



(6.4) 



lAK 



dD^ 



1/(01 

|5[(n,C)| 



|dC| ^ const^ / \C\^\dC\ ^ const(/)^(l + e)^ < 



const (/) 



dDe 



Due to assumption (j3.18p of Lemma 13.21 the integral I2 from (|3.17p is equivalent to zero. 
When estimating I^ and /4 from ()3.17p we fix some independent Eq > and integrate 
separately over Ds^^\D^ and C\D£q: 



|/:^l ^ 



D.o\De 



f'iC)eMtS{u,C)) 



S'^iuX) 



dReCdlmC + 



C\D,o 



/'(C)exp(t5(n,C)) 



S'^{u,C) 



dReCdImC ^ 



eo 

^const(/) f -^dp + const //" |/^(C)||C^|dReC(iImC ^ 



const (/) 



C\D,, 
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1/4 I ^ 



Deo\De 



fiC)expitS{u,0)S''{u,C) 



(S'AuXW 



'( 



+ 



dReCdImC+ 



C\D,g 



f{OeMtS{u,C))S'l^{u,0 



iS'AuX))' 



dReCdlmC ^ 



^ const (/) / ^dp + const It \f{Q)\\C,^\dRQ(,dlmQ ^ 



const (/) 



Setting finally e = ^^_^^^|^i/4 yields 



I{t,u) ^ 



const(/) 



uniformly on n G U. 

Case 2. u £ C\U and w from (|3.14|) satisfies ujq < j^ < 1 — uji for some fixed independent 
positive constants ujq and wi (i.e. the roots (q, C2 from (|3.14p are separated from T, defined by 
and the root C2 is separated from the origin) 

In this case the we can estimate 



s'/AuX) 



is'Au,oy 



constlCI^ , ^ r.^ 
^ ^ for C e dDp. 



Using these estimates and proceeding as in case 1, we obtain 

\Iint\ ^ const(/)e^, |/i| < const(/), I2 = 0, I/3I ^ const(/), I/4I ^ const(/)ln- 



Setting £ = jrm, we obtain that 



I{t, u) ^ const (/) 



ln(3+|t|) 



uniformly for the considered values of the parameter u. 

Case 3. u G C\U and j-^ < uq (i.e. the roots Co and C2 from ()3.14p lie in some neighborhood 
of T from (EH)) 



Lemma 6.1. For any t ^ to with some fixed to > and any cj > one of the following 
conditions holds 

(a) Q<u < ^ijijIiyTTi; 
(c) 3n: 



(l-f|t|)Tn + l/(2 + 27n+l) 



< o; < 



(i+|t|)'^"/(2+ 



2:^;;;^, where 7„+i = |7„ + i, 71 = |. 
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Proof. We note that 



7n+l , 1 

-, n — )• c«; 



2 + 27^+1 4' 

7n ^ 7n 



2 + 27„+i 2 + 27, 
71 ^1 



n 



2 + 272 8 

Thus the intervals from the cases (a), (b), (c) Vn G N cover the whole range < oj < +oo. D 

We will prove the result separately for three different cases depending on the value of pa- 
rameter OJ 

(a) < a; < 2e = j^TW^ 

In this case estimates (|6.2p . (|6.3p . (j6.4p hold and so the reasoning of the case 1 can be 

carried out to obtain that 

const(f) 

uniformly for the considered values of the parameter u satisfying 



» < - < jnw^- ("^^^ 



(b) OJ > ei/3 



(l + |t|)l/8 

In this case we estimate \Iint\ ^ const(/)e^. 

Further, we note that the derivative of the phase is estimated as 

15^^,01 ^^^^j^ for Ce5Z),. (6.6) 

Thus for Ii we obtain |/i| ^ const (/)-A3. 

In order to estimate the integral I^ we use the following estimate of the derivative 5^ for 
C G dDp when e ^ p ^ Sq: 



2 



S'AuX)\ > ^^ITTT^, if P <UJ, 



-CV-'s;i^ Id 



,3 



\S'{uX)\>'^^^, ifp> 



w. 



(6.7) 



It allows to derive j/sl ^ const(/)-^. 

Finally, we proceed to the study of the integral I4. We use the following estimates 



mu,c)y 



const I Cl^ 

P'^UJ'^ 



and 

J l/(C)K const(/)a;, if p < tj. 



(6.8) 
1/(01 ^const(/)/>, if p>L^. 



After integration we obtain the estimate I/4I ^ const (/) -273. 
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Setting finally e 



(l+|i|)3/8 



we obtain 



I{t,u) ^ 



const (/) 



(l + |t|)3/4 

uniformly for the considered values of the parameter u satisfying 



w > 



(6.9) 



(c) e'>'"+i < a; < 2e"'" , where e 

In -^ 1) 



(l + |t|)V8- 



We proceed similarly to the case (b). Evidently, lint can be estimated \Iint\ ^ const(/)e^+^". 
Employing the estimate ()6.6p we obtain |Ii| ^ const (/) i^"^-^ ■ 

Using ()6.7p in order to estimate I3 we obtain I/3I ^ const (/) ^7,|^^i . 
Finally, to estimate I4 we use ()6.8p and 



Q/l 



(S'^inX))' 



{ const I CI"* „ ^ , , 



to obtain I/4I ^ const (/)^ — !1L./£1 



Setting e 



(i+|t|)'/(^+2^"+ 



- yields 



\I{t,u)\^ 



const(/)ln(3 + |t|)) 



(l + |t|)3/4 

uniformly for the considered values of the parameter u satisfying 

1 2 



(1 + |i|)7n+l/(2+27„+i) 



<u; < 



1 + |t|)7n/(2+27„+i) ' 



.10) 



Finally, from Lemma l6. II it follows that we have proved the required estimate uniformly on 
the values of parameter u G C\U such that -^^ < wq. 

Case 4. n G C\U and j^ > 1 — wi (i.e. the roots C2, — C2 he in the wi-neighborhood of the 
origin) 

This case is treated similarly to the previous one. We denote ub = jr^. Then we use 
estimates 

|/(C)Kc(/)|^ + p|, 






constp , constpCb 
or |5^(n,C)| > |.|4 , 



(^;(^,c))^ 



const I Cl^ 
^ n or 



s'lci^X) 



iS'AuXW 



< 



const I Cl^ 



p^uj 



which hold for C, G Dp, to obtain the necessary estimates. 



D 



Proof of Lemma \4-l\ 1 • The proof of inequality ()4.10p repeats the proof of Lemma 13.21 The 
proof of inequality (j4.1ip also follows the scheme of the proof of Lemma 13.21 In this case 
we take Dg, to be the union of disks of the radius e with centers in the stationary points 
of 5(m, C) and in the point A. 
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For the case when A ^ T, where T is defined by (|6.ip . an estimate weaker than ()4.1ip 
can be obtained via a simphfied reasoning. Indeed, lint, as in (j3.16p . can be estimated 

\Iint\ ^ O {-(i^t\)^)- Using estimates ^^, ([Q]) and 

IC - A| ^ p for C G dDp (6.11) 



we obtain that | I^xt I ^ O ( (i,uy\i+s 3 ) • Setting e = ,,!".. ^^4 we get the estimate ()4.12p . 

2. In order to obtain estimates (j4.14p . ()4.15p we proceed according to the scheme outlined 
in Section [3j In this case the integral I2 does not annul. On the other hand, when the 
variable of integration belongs to T, the estimate (j4.1ip on the integrand of I2 is stronger 
than the estimate ()4.12p for the general case. Thus we obtain for ai(z,t) 

/ ^2 \ 



\h\ ^ o 




1 



1 



{l + \t\Ye 

Setting e = /-, , |!|-,i/4 yields the required estimate. The estimate ()4.15p is obtained simi- 
larly. 

3. We will give the scheme of the proof for estimate ()4.16p . The estimate ()4.17p is obtained 
similarly. 

We will prove ()4.16p by induction. Suppose that ()4.16p holds for all n = 1,2,..., A'". 
Then following the scheme of Section [3] and taking into account that d\{A^^ ■ /)(A) = 

(^^^^ • /)(A), we obtain forn = iV + 1: 

V(i + |t|)^+5r^iy \(i + |t|)^+5r^ie3^ 

1 




I/4I ^ o 



n+|t|)5+|r^l£3 



Setting e = ,,,-1/5 we obtain the required estimate. 

4. We represent Rz^ti^) as the sum of the following members 

Rz,tiX) = B{A + A^ + A^ + ...)■! + AB{A + A^ + A^ + ...)■ 1+ 

+ {A + A^ + A^ + .. .)AB{I + A + A^ + ...)■! = rI^^{\) + Rlt{\) + R:%tW. 

The convergence of the series at sufficiently large times follows from the estimate (j4.17p . 
Now let 

^l,(A) = ^^ + o(^),asA^oo. 

From (|i30D and (liTTD it follows that \qi{z,t)\ ^ y^^ijSS?? - From (liHD and (liTm we 
obtain that \q2{z,t)\ ^ f^.^'^zjl+i/r, - Finally, from (liT3|) . (liT6]) and (IJTTll it follows that 



(1+1*1)3/4+1/5 
(1+1*1)1/2+2/5 • 



1^3(2, t)| ^ i'i-i-mW2+2/5 • This yields the required estimate. D 
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